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We use the Brans-Dicke theory from the framework of General Relativity (Einstein frame), but 

now the total energy momentum tensor fulfills the following condition ^lx A "' M +T M "(0)j = 0. 

We take as a first model the flat FRW metric and with the law of variation for Hubble's parameter 
proposal by Berman pQ, we find solutions to the Einstein field equations by the cases: inflation 
(7 = —1), radiation (7 = ~), stiff matter (7 = 1). For the Inflation case the scalar field grows fast 
and depends strongly of the constant M 7= _i that appears in the solution, for the Radiation case, 
the scalar stop its expansion and then decrease perhaps due to the presence of the first particles. 
In the Stiff Matter case, the scalar field is decreasing so for a large time, <j> — > 0. In the same line 
of classical solutions, we find an exact solution to the Einstein field equations for the stiff matter 
(7 = 1) and flat universe, using the Hamilton-Jacobi scheme. 

PACS numbers: 04.60.Kz, 12.60.Jv, 98.80.Jk, 98.80.Qc 



I. INTRODUCTION 



The Scalar- Tensor theories have their origin in the 50' s. Pascual Jordan was intrigued by the appearance of a new 
scalar field in Kaluza-Klein theories, especially in its possible role as a generalized gravitational constant. In this way 
appear the Brans-Dicke theory, with the particularity that each one energy momentum tensor satisfy the covariant 
derivative [2J, T^ ul .„ = 0, where i corresponds to the i-th ingredient of matter content. Late years ago (1982), appear 
a new proposal by Barber [3J, known as self-creation cosmology (SCC) [H[S]. Since the original paper appeared in 
1982, more and more authors [MS] have worked the different versions of this theory in the classical fashion. By instant, 
Singh and Singh [lOj have studied Raychaudhary-type equations for perfect fluid in self-creation theory. Pimentel 
[TT] and Soleng [T21 [T3J have studied in detail the cosmological solutions of Barbers self-creation theories. Reddy 
[H] . Venkateswarlu and Reddy [IS], Shri and Singh QUE?], Mohanty et al. [18 , Pradhan and Vishwakarma [T9"ll2"U] . 
Sahu and Panigrahi [21], Venkateswarlu and Kumar [22] are some of the authors who have studied various aspects 
of different cosmological models in self-creation theory. These papers adapted the Brans Dicke theory to create mass 
out of the universes self contained scalar, gravitational and matter fields in simplest way. 

The gravitational theory must be a metric theory, because this is the easiest way to introduce the Equivalence 
Principle. However always is possible to put additional terms to the metric tensor, the most obvious proposal is a scalar 
field (j>. Recently Chirde and Rahate [33J investigated spatially homogeneous isotropic Friedman-Robertson- Walker 
cosmological model with bulk viscosity and zero-mass scalar field in the framework of Barber's second self-creation 
theory and found classical solutions, because is the simplest way to work this theory, because only take in account 
the energy-momentum tensor of usual matter and the scalar (f>. 

This work is arranged as follow. In section II we present the method used in general way, where we take the 
Brans-Dicke Lagrangian density and consider this theory in self-creation theory. In section III, employing the flat 
FRW cosmological model as a toy model with barotropic perfect fluid and cosmological constant and found classical 
solutions for three epoch in our universe, inflation phenomenon (7 = —1), radiation (7=3) and stiff matter (7 = 1). 
In Section IV we construct the Lagrangian and Hamiltonian densities for the cosmological model under consideration 
and are presented other class of classical solutions using the Hamilton-Jacobi approach. The section V is devoted to 
the conclusions of the work. 
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II. SELF CREATION COSMOLOGY IN GR 

The Lagrangian density in the Brans-Dicke theory is 



16tt 

where L mat ter = 167rp, so making the corresponding variation to the scalar field and the tensor metric, the field 
equations in this theory become 

fl-^r^,,. + ^nV = Q, (2) 
r<* _ i g ^R = 5E T <* + ^ - -g^W,*) + - - g^nV) , (3) 

where the left side is the Einstein tensor, the first term in the right side is the corresponding energy momentum tensor 
of material mass coupled with the scalar field <f> . The second and third term corresponds at energy momentum 
tensor of the scalar field coupled also to 4 > ~ 1 - Both equations are recombined given the relation 

U<t> = 4ttAT, (4) 

where A is a coupling constant to be determined from experiments. The measurements of the deflection of light 
restrict the value of coupling to |A| < 10 _1 . In the limit A — > 0, the Barber's second theory approaches the standard 
general relativity theory in every respect. Dip = 0:^ is the invariant D'Alembertian and T is the trace of the energy 
momentum tensor that describes all non gravitational and non scalar field matter and energy. Taking the trace of the 
equation Q and then substitute in Q, we obtain the following wave equation to <fi 

^=(3rk) T - < 5 > 

comparing both equations Q and (jHJ), we note that A = ^+ 2 u)) > wnere <^ is a coupling constant. 
Now the Einstein equation can be rewritten as 

Qa/3 _ S^rpc^M) + A T aW) _ T a/3(T)^ ^ 

where 

In Brans-Dicke theory, each one energy momentum tensor satisfy the covariant derivative, T M " - u — 0. In the self 
creation theory, we introduce that the total energy momentum tensor is which satisfy the covariant derivative, 
T ^(T) v = 0j namely; 



T^ T ) , p = 0, 



= 0, (7) 



which imply that [Q Q/3(T) ] = ^ [Q a/3(T) ], where Q Q ^ T ) = </>T Q/3 ( T ). This equation is the master equation which 
gives the name of self creation theory, because the covariant derivative of this tensor have a source of the same tensor 
multiply by a function of the scalar field <f>. 

III. FRW IN SELF CREATION THEORY 

We apply the formalism using the geometry of Friedmann-Robertson- Walker 

dr 2 



ds 2 = -N 2 (t)dt 2 + A 2 (t) 



1 — KT 2 



r 2 d<2 2 +r 2 sen 2 (6>)d^ 2 



(8) 
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where N is the lapse function, A is the scalar factor. Now we solve the equations (|3j), Q and Q, with the aim to 
find solutions to Density p(t), Scalar factor A(t) and the Scalar field 4>(t). 

Taking the transformation / = = 4^ and using T a/3 as a fluid perfect, first compute the classical field equations 
|3]), together with the barotropic equation of state P — 7/3, this equation become 

(A'\ 2 A' 4/ 1 (#\ 2 . p,„ k 
"{A +3 A6-2 UJ {6 -% +3 ^ = °' (9) 



A" A 12 A' ft 
A A 2 AS 



the equation Q become 



1 ^ p K 

+ +8 + _ =0 , 



(10) 



4ttA(1-37)- 



(11) 



The covariant equation ([7]) or conservation equation to the total energy momentum take the following form 

A ' / J.I \ 2 



3w- 



A' 



A' . p < 

+ 247T— (l+ 7 )^ +w - 



— 57T— — + 87T— = 0, 



(12) 



thus, the system equations to be solved are ([9^ - (12). 

We write the term (jt^j + -jt, m equation (9), using equations (10) and (11), and after some algebra we have 



A" A' <h' 

3— 3—— + w 

A A <b 



An 



A(l-3 7 ) + (l + 3 7 ) 



(13) 



Equation ( 12 ) can be rewritten as 



A" A' 6' 
-3— + 3—— - lj 
A A <h 



A' 
! A 



' o A' 

-w + 247r(l + 7 )^— +8tt 

1 A 



0. 



(14) 



and using the equations (11) and (13) in (14), then we have the master equation for solve the energy density of the 
model, as 



2tt [2(1 + 3 7 ) + A(l - 37) (3 + 2w)] 4^- + 24vr(l + 7) 4 4" + 87r 

(p m A 



= 0, 



defining the function F = 4, we have 



(15) 



dr 



Ln 



F A 3 ( 1 +7)0i[2(l+3 7 )+A(l-3 7 )(3+2c J )] 



who solution is 



p = M 7 A- 3 ( 1+ ^/, /? = [2 - Aa ] 



= 0, 



ag = 3 + 2w. 



(16) 



(17) 



This equation is equivalent to General Relativity expression 0] with the addition the last factor representing the self 
creation cosmology. We note that Ao is a new free parameter, however if Ao = A = 3+ 2 2cj , then (3 = so equation 
(17) became the usual solution to General Relativity. On the other hand for a photon gas 7 = 5, we have (3 = 
and equation (17) reduce to its General Relativity expression p = /9o(A/A ) -4 which is consistent, because in the 
radiation epoch there was not interaction between photon and the scalar field. 

Now we focus on finding the scalar field and scale factor 4>(t), A(t). Berman and Gomide fj] proposed the following 
law of variation for Hubble's parameter 



H = - = DA~ n , 
A 



(18) 
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where D and n are constants, 



dt ' 



using the following definition 



AA 



where q is the deceleration parameter. 
Using ( 18 1 into ( 19 ) we have 



1. 



q+1, 



(19) 



(20) 



the relation (20) imply that q as a constant. The sign of q indicated whether the model inflates or not. The positive 
sign of q i.e. q > correspond to standard decelerating model, whereas the negative sign — 1 < q < for < n < 1 
indicates inflation [24 . Many authors (see, Singh et al. 6, 8, 9J), have studied flat FRW and Bianchi type models by 



using the special law for Hubble parameter that yields constant value of deceleration parameter. The equation ( 18 ) 
writes as 



AA 



n-l 



D. 



Then solving the equation (21) we obtain the law for average scale factor A(t) as 

A(t) = Cie Dt , n = 0, 



(21) 
(22) 



A(t) = [nDt + C 2 ] 1/n , n^O, 



(23) 



where Ci and C 2 are constants of integration. Equation (23) implies that the condition for the expansion of the 
universe is n = q + 1 > 0. 

With the equation (17) and the Berman's law (23) we can find the solution to scalar field 4>(t), we can use C 2 = in 
our next calculations, we may recover the general results by the substitution nDt = C 2 + nDt'. Taking account again 
the classical field equations with A = and n = 0, adding equation ^ and ( 10 ) with the gauge TV = 1 we get 



.A 



A 



10 



.A, 
A< 



16tt(7-1)^ =0, 



(24) 



so let's solve the equation (24) for the following cases: 



Case I 



Inflation (7 = —1), then the equation (24) is written as 

2 



4+ 4 
A 



A 



A, 

•a; 



- 16tt- = 0, 



(25) 



now inserting (23) and (17) into (25), after some algebra we have 

"6-2n 



t z (j>+ -t<£ + 
n 



= leTrM-it 2 ^, 
i = ctV, 



t 4 (j) + &t(f> + b</> = ct z p , (26) 

where j3 = 2 — Ao(3 + 2us), the equation (26 1 has the form x 2 y" + axy' + by = cx n y m . We solve this equation for the 
case j3 = 1, so the equation (26) is rewritten as 



t>+ ~t(f> + 
n 



6 - 2n 



leTrM^t 2 ^ = 0, 



t 2 (f) + at</> + (ct 2 + b)4> = 0, 



which solution is [25] 



0(t) = t 1 ^ [CiI„(Vrt) + C 2 K„(^t)], 



5 



2 x 10 105 
1.8 x \O m ■ 
1.6 x 10 105 
1.4 x 10 lw 
1.2 x \O m ■ 
ij)(t) 10 x I0 m 
8 x 10 104 
6 x 10 1M 
4xl0 1M 
2 x 10 1M ■ 



/ 



/ 



/ 



/ 



/ 



10 x 10" 3S 3 x 10" 3i 5 x 10" !i 7 x 10" !i 10 x 10" !i 
f 



Figure 1: The shape of Scalar field is the same for < n < 1, the value of the constant M_i tells us how the scalar field grows, 
curve reaches the highest and lowest altitude for large and small values of the constant M_i respectively. In this case n = 0.8 
and M-i = 10 75 



wherel„, are a modified Bessel functions of first and second kind respectively, and a = — , b 



6-2n 



, c = 167rM_i, 

4b. Now we are interesting in the inflation time (10 _36 s ~ 10 _34 s), K„ diverges near the origin, 



then C2 = 0, so the solution to (271 is 



^(t)=t^[C 1 I„(^t)]. 



(27) 



Figure ([!]) shows the shape of the scalar field <fi in the inflation time. The behavior of the scalar field is the same for 
< n < 1 and depends strongly of the constant M 7= _i that appear in the solution (17 1 and (27). 
This solution was found with fi = 1, so in this case the physical quantities take a form 



P = p = CiM_it-«-[I I ,(Vct)]. 



(28) 



Case II 



Radiation (7 = ~). In this case the equation (24) is written as 



„A Ik 
A ( A 



.A4> 6 16 P 
>A6 + 6-T 7: 6 =0 > 



(29) 



16 M1/3 4 
t 2 + at0 + b0 = dt m , 



now ( 23 ) and 4 into ( 29 1 , after some algebra we have 

a- 5 • T6-2n 
n 

where a = b = S=£s d = ^r^, and m = 2s=i. 

n' ' 3 JnD] 4 / 11 1 n 

We need to solve the equation (30), so taking account the following transformation t = e x , then (30) is rewritten as 



(30) 



dx 2 1 dx 



de n 
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where ai = a — 1, now using the next transformation 



we 2 



So the equation (31 1 is now 



dx 2 



tu = de mx+ 5 aix , 



f= b- 



(32) 



(33) 



The solution of equation (33 1 depend the sing of the constant /, but remembering that < n < 1, so / 
then the solution is 

oj = Ci cosh(kix) + C2 sinh(kix) 



d 



c - /s sinh[ki(x - 



(^) 2 <o, 



(34) 



where xo is an arbitrary number and f = — kf < 0. Now solving the integral that appear in equation (34) with 
o" = m+ y, we have the solution to u is the following 



d 
kT 



to(x) — Ci cosh(kix) + C2 sinh(kix) + 

e — (ki-o-)xo 



! -k 2 2(k 1+ a) 



2(ki - a) 



a kix 



(35) 



Finally using the transformation t = e x , with x = In |t|, then the equation (32) became 



c/>(t) = CiCOsh(k! ln|t|)t" 



+ 



ki 



ki 



k? 



C 2 sinh(ki ln|t|)t-5 ai 

(<T+ki)x 



2(ki 



d 
"kT 



g— (ki— cr)x 

2(k x - a) 



(36) 



Equation ( 36 ) represent the solution to scalar field in the Radiation time and depend of the constants that appear 
there. Where 01 = $=r, a = § - d = ±^ J^L - 



3 " (nD) 4 / 



— , k\ = note that all constants depend on n and < n < 1. 



Figure ([2| shows the shape of the scalar field 4> in the inflation and radiation epoch, in the inflation time <fi 
is growing, but in the transition between inflation and radiation epoch, the scalar field slows its expansion and then 
decrease in the radiation time, perhaps due to the presence of the first particles. In the transition we must consider 
that the constant M 1 must change for each stage of the universe. 



Case III 



Stiff matter (7 = 1), then the equation (24) is written as 

2 



A 



A. 
'A, 



= o, 



(37) 



substituting equation ( 23 ) in ( 37 1 , yielding 



with a 



t 2 4> + -t</> + 



6-2n 



6 - 2n 



= 0, 



t 2 <p + &t<j) + b<p = 0, 



with the following solution 



cj>{t) = |i|^( Cl |t|" + c 2 \t\-»), (1 - a) 2 > 46, 



(38) 
(39) 



where c\, ci are constants and /i 



n-l 
2n ■ 



Then the solution to equation ( 38 ) is 



0(t) ^Clltl 1 "" +C 2 |t|-n. 



(40) 
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lntp(t) 
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Figure 2: Log(0)vst. The scalar field is increasing in the inflation epoch, but in the radiation time the scalar field decrease. 



Figure (|3| shows that the contribution of the scalar field at the beginning was significantly, but decreasing exponentially 
along the time, also we note that <fr decreases more quickly than in the radiation time. 
Finally, the physical quantities take the form: 



/9 = Mi(Dnt)T( Cl t 1 -5 +C2ir = ) 



(41) 



where /3 



Aq(3+2m) 
2 



On the other hand, remembering that our results are based on the Berman's law pQ, we obtain an expression for 
the average scalar factor (equations 22 23), taking the equation (23 ((case n ^ 0) to solve the corresponding set of 
equations. Figure Q shows the behavior of scalar factor to different values of n (0 < n < 1). Note that when n 
increases, the scalar factor grows rapidly, i.e. the slope of the curve is large, but for certain values of n, the slope 
of the curve approaches a plane curve. Remembering I? is a constant and in this case has the role of slope of the 
curve. With the observational data in your epoch, Berman [T] calculated this value for dust matter, yielding to the 
value D = 5.0 x io 28 ^S-. Considering the shape of the scale factor for some values in the n parameter, we say that 
the expansion scenary today corresponds for a special value to n, in the branch [0.3, 1) figure Qb)). These intervals 
can be modified when one calculate the value to the constant D considering the observational data today. This 
calculation will be part of forthcoming paper. 



IV. LAGRANGIAN AND HAMILTONIAN DENSITIES IN SCC 



In the previous section using the Berman's law, were solved the Einstein field equations. Now we will use the 
classical approach to find solutions to (p,A,cf>). The corresponding Lagrangian density using Q into ([I]) 

A0A 2 A 2 A<P uA 3 6 2 q 
C = 6^ + 6^ ~^- + 167rA 3 pN, (42) 
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Figure 3: The black and blue line represents the behavior of the scalar field in the stiff matter epoch. The red line represents 
the solution to the radiation epoch obtained in the previous section, we note that in the transition between radiation and stiff 
matter, the scalar field decreases more quickly than predicted radiation curve. For very large times 4> — > 0. 



the momenta are (IT g = ^) 

Acj)A A 2 <j) A 2 A wA 3 

A= 6(3 + 1)^ W + - An A), (43) 
^^^(An A -2 W , (44) 

when we write the canonical Lagrangian density ^canonical = n q q — NH, we obtain the corresponding Hamiltonian 
density as 

n = 4— 7[-6</» 2 n| + uA 2 u 2 A + 6A^n A n 

-l92wa A 6 ^p], a = (3 + 2w), (45) 



A. Classical scheme: Hamilton-Jacobi equation 

Using the gauge N = 12ao0A 3 , and the transformation in the momenta U q = ||, where S is known as the 
superpotential function, with this, the Hamiltonian density is written as (we include the equation (171 in this last 
equation) 



as 



-co A 2 



dS 
dA 



6Acj) 



dS 
dA 



-1927ra M 7 A" 3 ^- 1 V (/3+1) = 0, 



(46) 



This is the differential equation in the Hamilton-Jacobi theory, equation ( 46 ) , which can be solved for general 
barotropic fluid using the method of Lagrange-Charpit [2l)H2"8"] . solutions that will be presented elsewhere. So in the 
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Figure 4: (a) Behavior of the scalar factor in the inflation and radiation time, we note that for 0.6 < n < 1 is when the 
universe has a rapid expansion, (b) after radiation time the shape of the scalar field is modified so the expansion scenary today 
correspond to 0.3 < n < 1. D is the slope of the curve, but this constant have different values for each age of the universe. 
Berman p] calculate this value for dust matter (7 = 0) with D = 5.0 x 10 28£ ^-. 



follow, for simplicity we shall use the stiff matter case, 7 = 1 (so (3 + 1 = A °( 3 + 2 ") ). Using S(A, 
we have 



= S 1 (A) + S 2 (0), 



-60 



2 / dS 2 



! .2 fdSiV . , AdSA /dS 2 
to A 2 |^|+ 6A0 ' w " 



dA 



dA 



-1927ra Mi(/) (/3+1) = 0, 



(47) 



identifiying A^- = Ci = constant, then (471 can be solve as 



(48) 



using (|43|) and (|44|) we find that c\clq = 4 + §4, then (|48|) can be determined as 



0(t) = </>osech' 9 + 1 (— c -\/3aot) 



(49) 
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so, the solution for the scale factor A, become 

A(t) = A cosh/JTr (-coV3&ot) exp 



cp a t 

.0 + 1 



(50) 



V. CONCLUSIONS 



In this paper we have investigated FRW cosmological model of the universe in the framework of Barber's second 
self-creation theory, we obtain a solution for three cases: inflation (7 = —1), radiation (7 = |) and stiff matter 
(7 = 1), in the first case the scalar field was an increasing function and depends on the constant M_i, we note that 
in the transition between inflation and radiation epoch, the scalar field stop its expansion and then decreases. This 
behavior is due to the constant M 7 that appear in the solution must change for each stage of the universe and we 
must also consider the presence of the first particles in the radiation time, in the stiff matter time the scalar field 
decrease exponentially and for very long time <f> — > 0. Such that currently, the scalar field contribution is very small. 
In addition we found a new parameter Ao ^ A = 3+ 2 2 ^ more general than barber's parameter, but when Ao = A the 
solution of density (17) is the same that in general relativity, so there is no contribution of the scalar field. Our results 
are preliminary with this and will depend on how to adjust the constant with current observations. On the other 
hand, the classical solution found under the Hamilton- Jacobi approach (stiff matter) have an structure more general. 
Moreover, the feeling is that in some approximation, this class of solution could be tied with the Berman's law. 
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